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ABSTRACT 

A characterization  is  given  for  the  class  of  all  closed  loop  eigen- 
vector sets  which  can  be  obtained  with  a given  set  of  distinct  closed 
loop  eigenvalues  using  state  feedback.  It  is  shown,  furthermore,  that 
the  freedom  one  has  in  addition  to  specifying  the  closed  loop  eigen- 
values is  precisely  this:  to  choose  one  set  of  closed  loop  eigenvectors 
from  this  class.  Included  in  the  proof  of  this  result  is  an  algorithm 
for  computing  the  matrix  of  feedback  gains  which  gives  the  chosen  closed 
loop  eigenvalues  and  eigenvectors.  A design  scheme  based  on  these  results 
is  presented  which  gives  the  designer ^considerable  freedom  to  choose  the 
distribution  of  the  modes  among  the  output  conponents.  One  interesting 
feature  is  that  the  distribution  of  a mode  among  the  output  components 
can  be  varied  even  if  the  mode  is  not  controllable. 
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I.  INTRODUCTION 

One  of  the  most  popular  techniques  for  altering  the  response  charac- 
teristics of  a control  system  is  the  application  of  linear  state  variable 
feedback.  In  the  past  decade,  considerable  effort  has  been  made  to  under- 
stand exactly  what  feedback  has  to  offer  and  what  its  limitations  are. 

The  fact  that  one  can  use  state  feedback  to  assign  the  closed  loop 
system  any  desired  self  conjugate  set  of  eigenvalues,  provided  that  the 
open  loop  system  is  controllable,  is  a well  known  and  commonly  used  result. 
For  single  input  systems,  this  result  is  single  to  derive  and  has  been 
known  for  some  time.  Eigenvalue  placement  in  multi-input  systems  was 
studied  by  Lagenhop  [1],  Popov  [2],  Wonham  [3J,  Simon  and  Mitter  [4],  and 
Brunovsky  [5],  Nonham  was  the  first  to  prove  that  this  property  of  state 
feedback  also  applies  to  controllable  multi-input  systems. 

Eigenvalue  Assignment  in  multivariable  systems  is  still  not  well 
understood,  however.  Unlike  the  single  input  case,  specification  of  closed 
loop  eigenvalues  does  not  define  a unique  closed  loop  system.  Many  eigen- 
value placement  routines  offer  little  freedom  to  exploit  this  nonuniqueness. 
Exceptions  are  algorithms  which  allow  one  to  specify  a number  of  conponents 
of  the  closed  loop  eigenvectors  [6],  [7],  and  algorithms  designed  to  avoid 
large  feedback  gains  ([8]  chapter  6,  [9]). 

The  purpose  of  this  paper  is  to  identify  the  freedom  offered  by  state 
feedback  beyond  specification  of  the  closed  loop  eigenvalues  for  the  case 
in  which  the  desired  closed  loop  eigenvalues  are  distinct.  In  section  2 
of  the  paper,  a characterization  is  given  for  the  class  of  all  closed  loop 
eigenvector  sets  which  can  be  attained  with  a given  set  of  distinct  closed 
loop  eigenvalues.  It  is  shown,  furthermore,  that  the  freedom  one  has 
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beyond  specifying  the  closed  loop  eigenvalues  is  precisely  this:  to  choose 
one  set  from  the  class  of  "allowable"  sets  of  closed  loop  eigenvectors. 
Included  in  section  2 is  an  algorithm  for  computing  the  feedback  gain 
matrix  which  gives  a closed  loop  system  with  the  desired  eigenvalues  and 
eigenvectors.  Section  3 illustrates  one  use  of  eigenvector  assignment 
in  control  system  design.  Eigenvectors  may  be  selected  to  adjust  the 
distribution  of  the  modes  among  the  output  components . One  interesting 
feature  is  that  the  distribution  of  a mode  among  the  output  components 
can  be  varied  even  if  the  mode  is  not  controllable. 

Throughout  the  paper,  the  following  notation  will  be  used:  R,  C will 
denote,  respectively,  the  field  of  real  numbers  and  the  field  of  complex 
numbers.  For  a matrix  M which  has  n rows  of  numbers  in  R (C),Span  {M}  will 
be  the  subspace  of  R^C11)  spanned  by  the  columns  of  M,  while  Ker  {M} 
represents  the  Kernel  of  M (It  should  be  clear  from  the  context  whether 
the  operations  are  over  R or  C).  If  a property  is  said  to  hold  for  icn, 
this  means  that  it  holds  for  all  i satisfying  lgign.  For  a complex  number 
c,  c will  represent  its  complex  conjugate.  A set  is  described  as  self 
conjugate  if  the  complex  conjugate  of  each  of  its  menfcers  is  contained  in 


II.  SIMULTANEOUS  ASSIGNNCNT  OF  CLOSED  LOOP  EIGENVALUES  AND  EIGENVECTORS 
Consider  the  closed  loop  state  equation 
x(t)  • (A  ♦ BF)x(t)x(t) 

obtained  by  applying  linear  state  variable  feedback  u(t)  ■ Fx(t)  to  the 

system 

x(t)  ■ Ax(t)  + Bu(t) 

where  x(t)eRn#  u(t)€Rm,  rank  B « m . To  aid  in  the  development  of  the 
results,  we  associate  with  each  number  XeR  (C)  the  matrix 

Sx  - [XI-A  B] 

and  a compatibly  partitioned  matrix 


whose  colunns  constitute  a basis  for  Ker  {S^}.  It  is  quite  easy  to  show 

that  the  columns  of  N^  are  linearly  independent  if  B has  linearly  indepen- 

* 

dent  columns,  and  that  N^*  = N^  . 

It  is  well  known  that  F may  be  chosen  to  yield  any  self  conjugate  set 
of  closed  loop  eigenvalues  provided  that  the  (A,  B)  is  controllable;  that 
is,  provided  that 

rank  [B  AB  . . . An-1B]  - n 

or,  equivalently  (see  (8]),  that  rank  ■ n for  all  XcC.  Unless  m * 1, 
however,  specification  of  a set  of  closed  loop  eigenvalues  does  not  imiquely 
define  F.  It  is  a simple  matter  to  show  that  F is_  uniquely  defined,  if  it 
exists,  by  the  selection  of  a set  distinct  eigenvalues  together  with  a 
corresponding  set  of  eigenvectors.  Hence  the  design  freedom  available  beyond 


eigenvalue  selection  is  that  of  choosing  one  of  the  possible  closed  loop 
eigenvector  sets. 

Proposition  2.1,  the  main  result  of  this  paper,  treats  the  case  where 
the  closed  loop  eigenvalues  are  distinct.  For  this  case  necessary  and  suf- 
ficient conditions  for  the  existence  of  F which  yields  prescribed  eigen- 
values and  eigenvectors  are  given. * The  proof  includes  a procedure  for 
confuting  F. 

Proposition  2.1:  Let  {A^,  ien}  be  a self  conjugate  set  of  distinct  conplex 

numbers.  There  exists  a matrix  F of  real  numbers  such  that  (A  ♦ BF)  v^ 

* *ivi  and  0,1  ^ the  following  three  conditions  are  satisfied 

for  ien. 

1.  vectors  v^  are  linearly  independent  vectors  in  C° 

* * 

2.  * vj  whenever  A^  • X^ 

3.  v.  e Span{N  }. 

1 Ai 

If  F exists  and  rank  B * m,  then  F is  unique . 

Proof:  (Sufficiency) 

Suppose  that  v^  ien, are  chosen  to  satisfy  the  three  conditions  stated 

in  the  proposition.  Since  v.  e Span{N.  } for  ien  (condition  3),  then  v. 

1 i 1 

can  be  expressed  as  v.  » N..  k.  for  some  vector  k, eR™ CC111) , which  iiqplies 

X X X . 

that 

(X.  I - A)v4  ♦ BMj  k.  * 0. 

If  F is  chosen  so  that  -M^  k * Fv^,  then  [A^I  - (A*BF)]v^  = 0.  What 
remains  in  the  proof  is  to  show  that  a matrix  F of  real  numbers  satisfying 
F[Vj  v2  . . vn]  - [Wj  w2  . . wn]  (2.1) 

*i  ■ -\ki 

t Although  not  presented  in  this  form,  this  result  was  obtained  independently 
and  used  implicitly  by  H.  Kimura  in  his  treatment  of  pole  placement  using 
output  feedback  [10]. 
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can  always  be  constructed. 

If  all  n eigenvalues  are  real  numbers,  then  v^,  w^  are  vectors  of 
real  numbers  and  the  matrix  [Vj  v2  . . . vnJ  is  nonsingular.  For  this 
case 

F ■ [wj  w2  . . . w^lvj  v2  . . . vnJ_1. 

For  the  case  where  there  are  complex  eigenvalues,  assume  that  X^  ■ Xj, 

# 

The  second  condition  in  the  proposition  states  that  which  implies 

* 

that  Wj  ■ w2.  The  equation  which  must  be  solved  then  is 

F[ViR*jvn  vir-3vh  V]  - [w1R+jwn  w1R-iwn  (2*2) 
where  the  columns  of  V and  W are  v^,  i»3,  ...,n,  and  w^,  i»3,  ...,n, 
respectively.  Multiplication  of  both  sides  of  equation  (2.2)  from  the 
right  by  the  nonsingular  matrix 


\ -J*s 

0 

*1  +J*S 

0 

\ 

yields  the  equivalent  equation 

F*v1R  V1I  = 1W1R  W1I  W^* 

Clearly  since  v^,  icn,  are  independent,  the  columns  of  [v1R  Vjj  V]  are 
linearly  independent.  This  procedure  can  obviously  be  applied  for  all 
complex  pairs  of  eigenvalues. 

(Necessity) 

Necessity  of  the  first  two  conditions  follows  directly  from  elementary 
matrix  theory.  Furthermore  if  (A+BF)v^  * X^v^  then  [X^I  - A]v^  ■ BFv^; 


written  differently. 
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[X^-A  B] 


l 

■Fv. 


(2.3) 


Since  the  columns  of  K.  form  a basis  for  the  kernel  of  [X. I-A  B],  it 

Ai  1 

follows  that  v.  € Span{N.  }. 

1 l 

Since  A+BF  is  uniquely  defined  by  its  (distinct)  eigenvalues  and 
eigenvectors,  it  is  clear  that  F is  unique  whenever  B has  independent 
column  vectors. 

Q.E.O. 

Remark:  For  the  case  where  the  selected  eigenvalues  are  not  distinct, 

the  conditions  of  Proposition  are  sufficient,  and  the  last  two  conditions 
are  necessary. 

It  is  interesting  that  controllability  is  not  mentioned  in  this 

proposition.  The  equivalence  of  controllability  and  pole  assignability 

implies  that  the  three  conditions  in  Proposition  2.1  cannot  be  satisfied 

if  the  uncontrollable  eigenvalues  are  not  included  in  the  selected  set  of 

closed  loop  eigenvalues.  This  point  deserves  a few  words  of  clarification. 

If  (A,  B)  is  not  controllable,  then  there  exists  an  open  loop  eigen- 

T 

value  X.  and  a vector  p.  satisfying  p.S.  = 0.  This  implies  that  p.  is 

1 XX  Aj  X 

T 

an  eigenvector  of  A which  is  invariant  under  feedback;  i.e. 

(A+BF)TPi  * ATpA  * XJpi  . 

Now  for  any  X and  v e Span(N^)  there  exists  a matrix  F (see  the 
proof  of  Proposition  2.1)  satisfying  (A+BF)v  = Xv  . Hence 

Xj  (pjv)  * pT(A+BF)v  = X(pTv) 


tiimli'  MM  i M li  » 
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and  it  follows  that  « 0 for  all  \ i . 

The  important  point  is  this:  if  the  eigenvalues  of  the  uncontrollable 
subsystem  obtained  using  Kalman's  controllability  decomposition  (see  [11] 
or  almost  any  linear  systems  text)  are  distinct,  then  they  may  be  included 
in  the  set  of  selected  closed  loop  eigenvalues  without  violating  the  dis- 
tinct eigenvalue  assumption  of  Proposition  2.1.  Even  though  an  uncontrollable 
eigenvalue  is  invariant  under  feedback,  there  is  considerable  freedom, 

as  will  be  shown  in  section  3,  to  select  v.  e Span{N.  },  although  it  is 

1 Ai 

T 

clear  from  the  last  paragraph  that  v^  must  satisfy  p^  jf  0 . 


I 
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III.  APPLICATION  TO  CONTROL  SYSTEM  DESIGN 


While  the  overall  speed  of  response  of  the  closed  loop  system  is 
determined  by  its  eigenvalues,  the  "shape"  of  the  transient  response  depends 
to  a large  extent  on  the  closed  loop  eigenvectors.  In  the  following  para- 
graphs it  is  shown  that  through  eigenvector  selection  the  designer  has 
considerable  freedom  to  adjust  the  distribution  of  the  modes  among  the 
various  output  components.  In  this  sense,  he  may  "shape"  the  response 
characteristics  of  the  system. 

Consider  again  the  closed  loop  system 


x(t)  = (A+BF)x(t) 
together  with  an  output  equation 
y(t)  = Cx(t) 

where  y(t)eR  , rank  B=m,  rank  C=r,  m^r.  In  terms  of  the  (distinct)  closed  loop 
eigenvalues  A^,  ien,  and  eigenvectors  v^,  ien,  the  output  vector  is  given  by 

n T .A.t 

XCt)  = .^Cvi(pixo)e  1 

where 


[Pl  p2  . . . PnJT  = [Vj  v2  . . . vnJ_1. 
n A. 

In  other  words,  y(t)  = £ a.e  where  a.  is  proportioned  to  Cv. . Hence 

i=l  1 1 1 

T 

if  Cv,  = [1  0 ...  0)  then  the  ith  mode  appears  only  in  the  first 

T 

output  component.  If  Cv^  = [2  1 0 ...  0],  it  appears  in  the  first 

two  outputs  and  is  twice  as  large  in  y^t)  as  it  is  in  y2(t),  etc. 


Now  consider  the  freedom  available  to  assign  Cv  with,  vcN^,  where  A 

r r 

is  an  arbitrary  number  in  R (C).  The  columns  of  CN^  span  R (C  ),  meaning 
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that  cv  can  be  arbitrarily  assigned,  if  and  only  if 

* r + rank  [XI-A  Bj 


rank 


XI-A  B* 
C 0 


(3.1) 


To  see  this,  observe  that  there  exists  v s Span  such  that  Cv  = e if 
and  only  if  there  is  a solution  to  the  equation 

(3.2) 


T1 

HH 

1 

> 

03, 

1 — 

> 

1 

S 

1 

o 

J 

i — 

n 

,o 

L « J 

1 

A solution  to  this  equation  exists  for  every  vector  e iff 


fXI-A  B] 

XI-A  B O' 

rank 

i 

n 

i®_ 

= rank 

C 0 I_ 

The  rightmost  matrix  has  rank  r + rank  [XI-A  B],  If  (C,  A)  is  observable, 
this  can  be  stated  equivalently  as  follows:  Cv,  v c Span  N^,  can  be  arbi- 
trarily assigned  provided  that  X is  not  a zero  (as  defined  by  Rosenbrock 
[12])  of  the  transfer  function  matrix. 

Even  with  the  freedom  to  assign  Cv,  vcN^,  which  exists  if  equation 
3.1  holds,  complete  freedom  to  assign  Cv^,  ien,  where  are  closed  loop 
eigenvectors,  does  not  exist:  the  vectors  v^,  ien^,  must  satisfy  the  con- 
ditions of  proposition  2.1.  There  are  a few  obvious  constraints  on  the 
vectors  Cv^,  icru  It  is  clear,  for  example,  that  if  ra=r,  X is  not  an 
uncontrollable  eigenvalue,  and  (3.1)  holds,  then  Cv^=0  v^=0  which  is 

not  valid.  Considering  the  linear  independence  requirement,  it  is  also 
clear  that  the  number  of  modes  which  can  be  restricted  to  a single  output 
component  is  less  than  or  equal  to  n-r+1. 

The  case  in  which  an  eigenvalue  X^  is  uncontrollable  is  quite  interesting. 

In  the  last  section  it  was  shown  that  the  closed  loop  eigenvector  v^  must 
T T 

satisfy  p^v^  t 0,  where  A p^  = X^p^.  . If 


II 
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XI -A  B 0 


XI -A  B 


rank  C 0 0 * rank  C 0 

T T 

LpJ  oiJ  L Pi  o 


(3.3) 


then  there  exists  a vector  v^,  v^eN^  , such  that  Cv^  = 0 , p^v^  / 0 . 

This  implies  that  the  mode  can  be  eliminated  entirely  from  the  output 
response  (i.e.  made  unobservable)  provided  that  the  remaining  vectors 
Vj»  j€E»  0311  ke  c*105611  to  satisfy  the  conditions  of  Proposition  2.1. 

To  complete  this  section,  we  illustrate  one  design  scheme  which  allows 
a degree  of  transient  response  shaping.  It  is  assumed  that  the  uncon- 
trollable eigenvalues  are  distinct.  The  steps  of  the  procedure  are  as 
follows: 

1.  Select  distinct  closed  loop  eigenvalues  — the  set  must  include  un- 
controllable eigenvalues. 

2.  Select  e^,  the  desired  value  of  cv^,  for  each  eigenvalue.  This  choice 
should  be  made  on  the  basis  of  the  desired  distribution  of  modes  among 


output  components. 


3.  Compute  v^  satisfying 


[V-A  B][  vi 

L C oj  L w. 


If  X^  is  an  uncontrollable  eigenvalue,  the  solution  must  satisfy 
T 

PjV^  / 0.  This  inequality  may  be  forced  by  solving 


XiI-A  B 
C 0 


e.  , k 4 0 

l ' 


if  the  leftmost  matrix  has  rank  r+l+rank[X^I-A  BJ  . 
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4.  If  the  vectors  v^,  ien  are  not  linearly  independent,  alter  one  or  more 

of  the  vectors  e^,  icn_,  and  return  to  step  3;  otherwise  proceed  to  step  5. 

5.  Compute  the  feedback  matrix  and  the  closed  loop  system  matTix.  If  the 
transient  response  characteristics  are  not  satisfactory  alter  one  or 
more  of  the  vectors  e^,  according  to  the  nature  of  the  response,  and 
return  to  step  3. 


Example:  Consider  the  open  loop  system  with  matrices 

1-1.25  0.75  -0.751 

A = | 1 -1.5  -0.75 

1 -1  -1.25 


B « 


C « 


1 

0 

0 

1 

0 


0 

1 

1 

0 

0 


0 

1 


which  has  controllable  eigenvalues  at  -1.25,  -2.25,  and  an  uncontrollable 
eigenvalue  at  -0.5. 

To  illustrate  the  response  shaping  feature  of  the  procedure  suggested 

in  this  section,  it  is  assumed  that  the  objective  is  to  shift  the  controllable 

eigenvalues  to  -5.,  -6.,  and  to  obtain  a "reasonable"  response  for  the 

T 

initial  condition  [001],  which  represents  a disturbance  in  Xj(t). 

* 

Table  1 gives  three  different  closed  loop  systems  which  have  eigenvalues 
Xj  = -.5,  = -5,  = -6,  but  which  differ  in  their  eigenvectors. 


The  first  system  in  Table  1 was  obtained  using  a standard  modal  con- 
troller (see  [8],  chapters  5,6)  of  the  form 


* Computed  answers  have  been  rounded  to  three  significant  figures.  Numbers 
smaller  in  magnitude  than  lO-^  are  shown  as  zero. 


■ - 
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u(t)  * (“jgjpj  ♦ “282p2)x(t) 
where  the  following  equations  hold: 

j 

ATPj  * -1.2Sp1#  pjpj  « 1 

(A  + Bg1Pj)p2  = -2.25p2,  p2p2  = 1 . 

The  constants  c^,  o2  and  the  vectors  gj,  g^R11  were  chosen  to  shift  one 

eigenvalue  from  -1.25  to  -5.0,  and  then  the  second  eigenvalue  from  -2.25 

to  -6.0,  and  to  minimize  the  gain  required  for  each  shift.  Figure  1 shows 

the  response  of  this  system  to  a unit  disturbance  in  x^Ct) . 

The  response  of  the  first  system  is  not  satisfactory  because  of  the 

rather  large  slow  mode  component  in  y^t).  This  mode  is  uncontrollable, satisfies 

equation  (3.3),  and  is  eliminated  in  the  output  of  system  2 (Figure  2).  In  this 

system  e2,  e^  were  chosen  to  be  equal  to  cv*,  cv^,  where  v2,  v*  are  the 

eigenvectors  corresponding  to  A 2,  A^  in  system  1.  Note  that  the  entries 
1 

of  cv2  have  opposite  signs. 

To  further  illustrate  the  spirit  of  the  procedure,  it  is  assumed  that 
a positive  going  transient  in  y^(t)  is  desired  for  a positive  disturbance 

■ 

in  a x3(t).  This  is  accomplished  in  system  3 by  simply  choosing  e2  to  have 

T 

components  of  the  same  sign,  in  this  case  e2  = [.5  .5]  . Figure  3 gives 

the  response  of  the  third  system. 


J 


IV.  CONCLUSIONS 


Necessaiy  and  sufficient  conditions  for  simultaneous  eigenvalue,  eigen- 
vector assignment  have  been  given  for  the  case  where  the  desired  eigenvalues 
are  distinct.  The  corresponding  design  procedure  based  on  this  result  gives 
the  designer  considerable  freedom  to  select  the  distribution  of  modes  (whether 
controllable  or  not)  among  the  output  components. 

The  shortcoming  of  the  design  procedure  is  that  it  is  ad  hoc  in  nature. 
In  future  research,  an  attempt  will  be  made  to  develop  a systematic  design 
procedure  allowing  the  designer  to  "shape"  in  some  sense  the  average  tran- 
sient response  characteristics  of  the  system. 
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